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Anomalous properties in 
ferroelectrics induced 
by atomic ordering 

A. M. George, Jorge iniguez & L Bellaiche 

Physics Department, University of Arkansas, Fayetteville, Arkansas 72701, USA 



Complex insulating perovskite alloys are of considerable techno- 
logical interest because of their large dielectric and piezoelectric 
responses. Examples of such alloys include (Ba^rJTiOa, which 
has emerged as a leading candidate dielectric material for the 
memory-cell capacitors in dynamic random access memories 1 ; 
and PMZr^TiJOs (PZT), which is widely used in transducers 
and actuators 2 . The rich variety of structural phases that these 
alloys can exhibit, and the challenge of relating their anomalous 
properties to the microscopic structure, make them attractive 
from a fundamental point of view. Theoretical investigations of 
modifications to the atomic ordering of these alloys suggest the 
existence of further unexpected structural properties 3 and hold 
promise for the development of new functional materials with 
improved electromechanical properties. Here we report ab initio 
calculations that show that a certain class of atomic rearrange- 
ment should lead simultaneously to large electromechanical 
responses and to unusual structural phases in a given class of 
perovskite alloys. Our simulations also reveal the microscopic 
mechanism responsible for these anomalies. 

The class of Pb(Mg 1/3 Nb^)0 3 -PbTi0 3 (PMN-PT) and 
Pb(Zn 1/3 Nb M )0 3 --PbTi0 3 (PZN-PT) perovskite ferroelectric alloys 
have been reported 4 to show remarkably large piezoelectric 
constants around 2,000 pCN"\ These materials thus promise 
improvements in the resolution and range of ultrasonic and sonar 
listening devices 5 . Perovskite A(B'B"...)0 3 alloys are also of great 
interest, as demonstrated by the discovery of an unexpected 
monoclinic phase 6 in Pb(Zri_ x Ti x )0 3 . These two findings have led 
to a search for materials with even larger electromechanical 
response and/or currently unobserved phases. In particular, if a 
mechanism can be found that occurs in a large number of 
ferroelectric alloys, markedly enhances the piezoelectricity and 
the dielectric responses, and leads to unexpected structural 
features, that mechanism is likely to have large technological and 
fundamental implications. Here we report that such a mechanism 
exists and simply consists in rearranging in a certain way the atoms 
in a heterovalent alloy (that is, an alloy made of atoms belonging 
to different columns of the periodic table). More precisely,, we 
predict that materials made of the sequences PKSc^ F NbJJ., J(y 
Pb(Sc£r^^ alongyjthe 
[001] direction should show very large electromechanical responses, 
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and unusual monoclinic and orthorhombic phases for some values 
of the v parameter. The analysis of our results points to the different 
valence of the B atoms (Sc and Nb) as the main reason for the 
existence of these anomalous properties. 

Here, we use the numerical scheme proposed in refe 7 and 8, 
which consists of constructing an effective hamiltonian for the alloy 
from first-principles calculations, to predict the properties of 
PbCSc^NbJOj (PSN) structures. This effective hamiltonian 7 - 8 
contains a local-mode self-energy (expanded in even-order terms 
up to fourth order of the phonon local soft modes), a long-range 
dipole-dipole interaction, a short-range interaction between soft 
modes (quadratic order of the soft modes), an elastic energy 
(expanded to second order of the strain variables), and an inter- 
action between the local modes and local strain (second order of the 
soft modes and first order in strain). It also contains an energy term 
describing the effects of the atomic configuration on the phonon 
local soft modes (which are directly related to the electrical 
polarization) which can be written as 



AE 



(1) 



where i runs over all the cells and ; runs over the three nearest- 
neighbour shells of cell i. Here u, is the local soft mode in cell i, Z* is 
the Born effective charge associated with the local soft modes, and e^ 
is the unit vector joining the B site; to the B site i. The variables {o)} 
characterize the atomic configuration: a } = +1 ( - 1) indicates that 
there is an Nb (Sc) ion in cell;'. Sji is an alloy-related parameter that 
only depends on the distance between the B sites i and ;, and 
is derived from first-principles calculations 9 " 12 on small cells. 
Equation (1) indicates that AE can be viewed as the interaction 
energy between the dipole moment Z*u, associated with the site / 
and an internal electric field E ; - - afify induced by the B ions of 
sites ; on the site i. Moreover, we numerically find that the S ;7 
parameters all have negative signs. As a result, the radial electric 
field - afytgi acting on site i and induced by a Sc 3+ (Nb 5+ ) ion 
sitting on site ; is directed from the site i to the site ;" (respectively, 
from site ; to site i). This is consistent with an electrostatic picture of 
PSN since Sc (Nb) ions are negatively (positively) charged with 
respect to the average B-ion valence of 4+. The effective hamiltonian 
approach yields good agreement with direct first-principles results 
in PZTand PSN alloys 7 ' 8 . Previous works (see refs 7-8) have found 
that a linear rescaling of the simulation temperature often leads to 
good agreement with experiment. Adopting this approach again 
here, our temperatures are rescaled down by a factor of 2.5 so that 
the theoretical Curie temperature in disordered PSN is forced to 
match the experimental value 13 . The need for such rescaling may be 
due to higher perturbative terms or the rotation of the oxygen 
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Figure 1 Schematic representation of the ferroelectric structures considered here. The 
rnaterialsarePb(Sc^ 

Pb(Sc^Nb^)0 3 structures, Trie internal electric fields acting on the four B planes are 
indicated by arrows. 



octahedra, neglected in our effective model of PSN. 

Here , we use the total energy of our alloy effective hamiltonian in 
Monte Carlo simulations to calculate finite-temperature properties 
of some selected Pb(Sc,_ x Nb Jt )0 3 structures. We consider structures 
with the following sequence of four B-planes along the [001] 
direction (Fig. 1): a niobium-poor plane for which x = 0.5 - v, a 
second plane made of 50% scandium and 50% niobium (x = 0.5), a 
niobium-rich plane with x = 0.5 4- v , and a fourth plane similar to 
the second one. The B atoms are randomly distributed within each 
of these four planes. The studied structures only differ in the value of 
the parameter p, which is allowed to vary from zero (in the case of 
the disordered PSN alloy) to 0.5 (in the case in which the first plane 
is entirely made of Sc atoms and the third plane is fully occupied by 
Nb atoms). We use 12 X 12 X 12 supercells, implying that the 
sequence of the four different B(001) planes is, repeated three 
times, to get well converged results 7,8 . 

Figure 2a shows the cartesian coordinates (w*, viy and u 2 )— along 




Figure 2 Properties of the studied structures as a function of the v parameter at 20 K. 
a, Average cartesian coordinates u x , u r and u t of the local mode; b, the fa piezoelectric 
coefficient; c, the X33 dielectric susceptibility. Trie filled symbols in a represent the 
cartesian coordinates ^(LDA), u,(LDA) and u 7 (LDA) of the local modes predicted by 
direct first-principles calculations made using the local density approximation 9 " 12 , a.u., 
arbitrary units. 
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the [100], [010] and [001] directions, respectively — of the supercell 
average of the local mode vectors as a function of this parameter v at 
20 K. The structure for which v is null has u x = u y = u z ^ 0. This 
characterizes a ferroelectric rhombohedral structure in which the 
polarization is directed along the pseudo-cubic [111] direction, 
consistent with experiments on disordered PSN samples 13 . Inter- 
estingly, increasing v results in a strong decrease of u z whereas u x and 
tiy slowly increase and remain equal to each other This behaviour 
corresponds to a ferroelectric phase of monoclinic symmetry for 
which the polarization lies between the [111] and [110] directions 
(the one denoted M B in ref. 14). When v is larger than 0.44, u ? 
becomes null and u x and Uy reach their maximum value, indicating 
that the resulting phase is now of orthorhombic symmetry with a 
polarization lying along the [110] direction. Increasing the para- 
meter v thus leads to three different ground states and to a 
continuous rotation of the electrical polarization. To our knowl- 
edge, the predicted monoclinic M B ground state has never been 
observed in any perovskite material without the aid of external stress 
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Rgure 3 Properties of the structure associated with v = 0.375 as a function of 
temperature. Panels a, b and c display the same property as the corresponding panels of 
Fig. 2. Dashed iine in a shows one of the three cartesian coordinates (u x ) of the local mode 
in disordered (dis) PSN (the two other cartesian coordinates are nearly identical to the 
one displayed, and are omitted for clarity). 



or electric field. Similarly, we are not aware of any studies pointing 
to an orthorhombic ground state in PSN. The existence of these 
monoclinic and orthorhombic phases is due to the energy term 
described in equation (1). Figure 2a also shows the predictions of 
local density approximation 9 " 11 calculations performed on 20-atom 
supercells mimicking the structures corresponding to v = 0.25, 
v = 0.36 and v = 0.50. In these supercells, each B(001) plane is 
represented by a virtual atom corresponding to the composition in 
this plane 12,15 . The direct first-principles calculations agree with 
the predictions of our effective hamiltonian approach, thus 
further supporting the existence of these unusual monoclinic and 
orthorhombic phases. 

Figure 2b summarizes the effects of the parameter v on piezo- 
electricity and Fig. 2c the effects on dielectric response at 20 K. We 
found that the electromechanical coefficients most affected by the 
atomic ordering are the shear d M piezoelectric coefficient and the X33 
dielectric susceptibility, when expressing both piezoelectric and 
dielectric tensors in the orthonormal basis formed by a, = [100], 
a 2 = [010] and a 3 = [001]. More precisely, Fig. 2b shows that d u 
reaches a peak and remains at large values for a broad range of v 
centred around the monoclinic-to-orthorhombic phase transition. 
Figure 2c demonstrates that the studied atomic ordering simulta- 
neously results in a large dielectric response, as X33 achieves values 
above 1,000 at 20 K when v is greater than 0.3. 

We now investigate the finite-temperature properties of a struc- 
ture showing one of the largest electromechanical responses at 20 K. 
More precisely, we focus on the structure for which v = 0.375. 
Figure 3a shows the cartesian coordinates of the supercell average of 
the local mode vectors as a function of temperature for this 
structure, and compares them with those of the disordered PSN 
material. Each coordinate of the local mode in each structure is close 
to zero at high temperature, characterizing a paraelectric phase. At a 
temperature close to 373 K, the disordered material undergoes a 
transition from a paraelectric cubic phase to a ferroelectric rhom- 
bohedral structure (for which u x = u y = « z ), consistent with 
experiments 13 . The disordered PSN material then remains in the 
rhombohedral phase for lower temperature. On the other hand, the 
modulated structure with v = 0.375 adopts three different phases: a 
paraelectric tetragonal phase induced by atomic ordering (for which 
u x = u y = u z = 0) at high temperature, an orthorhombic ferro- 
electric phase (u x = u y =t0 and u z = 0) for temperature between 
400 K and 40 K, and the monoclinic ferroelectric M B phase 
(u x — u y > u z , with u z # 0) for temperature lower than 40 K. As 
shown in Fig. 3b and c, the existence of the orthorhombic-to- 
monoclinic phase transition results in huge electromechanical 
responses peaking around this transition and occurring over a 
broad range of temperature. In fact, d M and X33 are greater than 
1,500 pCN" 1 and 3,000, respectively, for any temperature lower 
than 100 K. 

We also numerically found (results not shown here) that con- 
tinuously increasing v from 0 to 0.44 leads to a continuous decrease 
of the orthorhombic-to-monoclinic transition temperature from 
373 K to 0 K. As a result, the temperature at which d M and X33 both 
peak depends on the value of the v parameter. This dependency 
could lead to the development of devices with electromechanical 
performances optimized for any temperature between 373 K and 
OK. 

The intriguing results of Figs 2 and 3 can be understood by means 
of equation ( 1 ). The internal electric field E ; - o } S^ acting on each 
B(001) plane is represented in Fig. 1. The atoms in plane 2 feel an 
internal electric field oriented along the [001 ] direction whereas the 
atoms in plane 4 feel an internal electric field along the [001] 
direction. It is straightforward to demonstrate that the magnitude 
of these fields is linearly dependent on the parameter v. (These 
features can be understood qualitatively by simple short-range 
electrostatic considerations: the difference of valence between Sc 
and Nb atoms leads to charged B(001) planes with ionic charges 
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linearly dependent on v. When electrostatic interactions up to the 
third-neighbour shells are included, these charged planes generate 
internal electric fields whose directions are those of Fig. 1 and whose 
magnitudes are linearly dependent on v.) Increasing v thus results in 
stronger internal fields with opposite directions. These strong 
opposite fields tend to suppress the polarization's component 
along the direction of compositional modulation. As a result, u z 
becomes smaller than u x and The studied structures with 
intermediate values of v thus first undergo a paraelectric-to- 
orthorhombic ferroelectric transition at high temperature before 
adopting the monoclinic M B phase for ground state. The modulated 
structures with the largest values of v have the strongest internal 
electric fields, which annihilate u z at any temperature. Conse- 
quently, such structures never reach the monoclinic phase and 
instead adopt the orthorhombic phase for ground state. The large 
values of d M and X33 shown in Figs 2 and 3 simply reflect the 
considerable change of u z when some parameters are slightly 
modified (v at a fixed temperature, or T for a fixed value of p), 
especially for structures at the borderline between the monoclinic 
and orthorhombic phases. In other words, the large electromecha- 
nical responses are consistent with the ease of rotating the 
polarization 716 . We have shown here results for the thinnest possible 
structures. For larger structures, with thicker layers, the 'unusual' 
electromechanical responses will be smaller due to smaller internal 
electric fields. 

On the basis of the electrostatic considerations discussed above, 
we expect that any alloy made of heterovalent atoms and with a 
rhombohedral ground state in its disordered form should have the 
structural, piezoelectric and dielectric anomalies displayed in Figs 2 
and 3, when the atomic ordering along the [001] direction is 
adjusted in a certain way. The atomically ordered structures dis- 
cussed here could be grown by means of a pulse laser deposition 
technique 17 or by using molecular beam epitaxy. □ 
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